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Abstract. We have used the radiation held from a star- 
burst population synthesis model appropriate for the 
brightest H n region of I Zw 18 to perform a photoion- 
ization model analysis of this object. We have investi- 
gated whether, with the observed nebular density distri- 
bution as revealed by the HST images and a total stel- 
lar radiation compatible with the observed UV flux, one 
could reproduce the constraints represented by the ob- 
served ionization structure, the He II A4686/H/3 ratio, the 
Ha flux and the electron temperature indicated by the 
[O m] A4363/5007 ratio. 

We have found that, even taking into account strong 
deviations from the adopted spectral energy distribution 
of the ionizing radiation and the effect of additional X- 
rays, the photoionization models yield too low a [O III] 
A4363/5007 ratio by about 30%. This discrepancy is sig- 
nificant and poses an interesting problem, which cannot 
be solved by expected inaccuracies in the atomic data. 
The missing energy may be of the same order of magni- 
tude as the one provided by the stellar photons or lower, 
depending on the way it acts on the [O in] A4363 line. 

Elemental abundance determinations in I Zw 18 are 
affected by this problem. Before the nature of the miss- 
ing heating source and its interactions with the nebular 
gas are better understood it is, however, not possible to 
estimate the typical uncertainties by which standard em- 
pirical methods are affected. 

Several side-products of our photoionization analysis 
of I Zw 18 are the following: 

— We have been able to reproduce the intensity of the 
nebular He n A4686 using a stellar radiation field 
consistent with the observed Wolf-Rayet features in 
I Zw 18. 

— We have shown that the bright NW H n region in 
I Zw 18 does not absorb all the ionizing photons from 
the central star cluster, and that about half of them 
are available to ionize an extended diffuse envelope. 

— The [O i] emission can easily be accounted for by con- 
densations or intervening filaments on the line of sight. 
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— The intrinsic Ha/H/3 ratio is enhanced by collisional 
excitation to a value of 3. Previous estimates of the 
reddening of I Zw 18 are therefore slightly overesti- 
mated. 

— We give ionization correction factors appropriate for 
I Zw 18 that can be used for more accurate abundance 
determinations in this object once the electron tem- 
perature problem is better understood. 

Key words: Galaxies: abundances - Galaxies: ISM - 
Galaxies: starburst - Galaxies: individual: I Zw 18 - Stars: 
early-type - Stars: Wolf-Rayet. 



1. Introduction 

The blue compact dwarf emission-line galaxy I Zw 18 is fa- 
mous for being the most metal poor galaxy known so far. 
Its oxygen abundance is about 2% the solar value, as first 
shown by Searle and Sargent (1972), and then confirmed 
by many studies (e.g. Lequeux et al. 1979, French 1980, 
Kinman & Davidson 1981, Pagel et al. 1992, Legrand et 
al. 1997, Izotov et al. 1997b, Vilchez & Iglesias-Paramo 
1998). Because of this, I Zw 18 has played an essential role 
in the determination of the primordial helium mass frac- 
tion. Also, due to its extreme properties, I Zw 18 has been 
a choice target for studies of star formation history in 
blue compact galaxies (Dufour & Hester 1990, Hunter & 
Thronson 1995, Dufour et al. 1996, De Mello et al. 1998, 
Aloisi et al. 1999), of the elemental enrichment in dwarf 
galaxies (Kunth & Sargent 1986, Kunth et al. 1995) and of 
the interplay between star formation and the interstellar 
medium (Martin 1996, van Zee et al. 1998). An important 
clue is the distribution of the oxygen abundance inside 
the H II regions (Skillman & Kennicutt 1993, Vilchez & 
Iglesias-Paramo 1998) and in the neutral gas (Kunth et 
al. 1994, Pettini & Lipman 1995, van Zee et al. 1998). 
Another clue is the carbon and nitrogen abundance (Gar- 
nett et al. 1997, Vilchez & Iglesias-Paramo 1998, Izotov 
& Thuan 1999). On the whole, there is a general consent 
about an intense and recent burst of star formation in 
I Zw 18 - which provides the ionizing photons - following 
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previous star formation episodes. How exactly has the gas 
been enriched with metals during the course of the evolu- 
tion of I Zw 18 remains to be better understood. 

Much of our understanding (or speculations) on the 
chemical evolution of IZwl8 (and other galaxies in gen- 
eral) relies on the confidence placed in the chemical abun- 
dances derived from the lines emitted in the H II regions. 
These are generally obtained using standard, empirical, 
methods which have been worked out years ago, and rely 
on the theory of line emission in photoionized gases. Pho- 
toionization models are, most of the time, used merely as 
a guide to evaluate the temperature of the low excitation 
regions once the characteristic temperature of the high 
excitation zones has been obtained through the [O III] 
A4363/5007 ratio. They also serve to provide formulae 
for estimating the correction factors for the unseen ionic 
species of a given element. 

Direct fitting of the observed emission line spectrum 
by tailored photoionization models provides more accu- 
rate abundances only if all the relevant line ratios are per- 
fectly reproduced by the model (which is rarely the case 
in model fitting history) and if the model reproducing all 
the observational constraints is unique. 

One virtue of model fitting, though, is that it permits 
to check whether the assumptions used in abundance de- 
terminations are correct for a given object. For example, 
there is the long standing debate whether so-called "elec- 
tron temperature fluctuations" (see e.g. Mathis 1995, Pc- 
imbcrt 1996, Stasinska 1998) are present in H n regions 
to a sufficient level so as to significantly affect elemental 
abundance determinations. If a photoionization model is 
not able to reproduce all the temperature sensitive line 
ratios, the energy balance is not well understood, and one 
may question the validity of abundance determinations. 
Also, photoionization models are a potential tool (see 
e.g. Esteban et al. 1993, Garcia- Vargas 1996, Stasinska 
& Schaerer 1997, Crowther et al. 1999) to uncover the 
spectral distribution of the ionizing radiation field, thus 
providing information on the ionizing stars, their evolu- 
tionary status and the structure of their atmospheres. 

These two points are a strong motivation for a pho- 
toionization model analysis of I Zw 18. There have already 
been a few such attempts in the past (Dufour et al. 1988, 
Campbell 1990, Stevenson et al. 1993). None of those mod- 
els were, however, able to reproduce the He II A4686 line, 
known to exist in I Zw 18 since the work of French (1980). 
The reason is that, in those models, the spectral distri- 
bution of the ionizing radiation was that of a single star 
whose radiation field was interpolated from a grid of plane- 
parallel, LTE model atmospheres for massive stars. Re- 
cently, Wolf-Rayet stars have been identified in IZwl8 
through the characteristic bump they produce at 4650 A 
(Izotov et al. 1997a, Legrand et al. 1997). Spherically ex- 
panding non-LTE model atmospheres for hot Wolf-Rayet 
stars with sufficiently low wind densities (Schmutz et al. 
1992) do predict an output of radiation above the He II 



ionization edge, which might, at least qualitatively, pro- 
vide a natural explanation for the narrow He n A4686 line 
observed in I Zw 18. Schaerer (1996) has, for the first time, 
synthesized the broad (stellar) and narrow (nebular) He n 
features in young starbursts using the Geneva stellar evo- 
lution tracks and appropriate stellar model atmospheres. 
He then extended his computations to the metallicity of 
IZwl8 (De Mello et al. 1998). 

In this paper, we use the emergent radiation field from 
the synthetic starburst model presented in De Mello et 
al. (1998) to construct photoionization models of IZwl8. 
One of the objectives is to see whether this more realistic 
ionizing radiation field permits, at the same time, to solve 
the electron temperature problem encountered in previ- 
ous studies. Former photoionization models predicted too 
low a [O in] A4363/5007 ratio, unless specific geometries 
were adopted (Dufour et al. 1988, Campbell 1990), which 
later turned out to be incompatible with Hubble Space 
Telescope (HST) images. The synthetic starburst model 
we use is based on spherically expanding non-LTE stel- 
lar atmosphere models for main sequence stars (Schaerer 
& de Koter 1997) and for Wolf-Rayet stars (Schmutz et 
al. 1992). These models have a greater heating power than 
the LTE model atmospheres of same effective temperature 
(see Fig. 3; also Schaerer & de Koter) 

The progression of the paper is as follows. In Section 2, 
we discuss in more detail the photoionization models pro- 
posed previously for I Zw 18 and show in what respect they 
are not consistent with recent observations. In Section 3, 
we present our own model fitting methodology, including 
a description of the computational tools. In Section 4, we 
describe the models we have built for I Zw 18, and discuss 
the effects of the assumptions involved in the computa- 
tions. Our main results are summarized in Section 5. 



2. Previous photoionization models of I Zw 18 

The first attempt to produce a photoionization model for 
IZw 18 is that of Dufour et al. (1988). Their observational 
constraints were provided by spectra obtained in an aper- 
ture of 2.5" x 6"of the NW region combined with IUE ob- 
servations yielding essentially the C in] A1909 line. Using 
Shiclds's photoionization code NEBULA, they modelled 
the object as an ionization bounded sphere of constant 
density n = 100 cm -3 and adjustable volume filling fac- 
tor e so as to reproduce the observed [O in] A5007/[O n] 
A3727 ratio. The ionizing radiation was provided by a cen- 
tral source of radiation, represented by the LTE model at- 
mospheres of Hummer & Mihalas (1970), modified to take 
into account the low metallicity of I Zw 18. Discarding the 
He II problem, they obtained a model that was reason- 
ably successful except that it had an ++ temperature, 
T(0 ++ ), marginally smaller than observed (17200 K) 
compared to the value of 18100 (+1100, -1000) K derived 
directly from their observed [O ill] A4363/5007 = (2.79 ± 
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0.35) x 10~ 2 (the errors quoted being 2a) []. This model 
was obtained for an effective temperature of 45000 K. 
These authors showed that, because of the dominant role 
played by Lya cooling in IZwl8, it was impossible, for 
the adopted geometry, to produce a model with notice- 
ably higher T(0 ++ ), by varying the free parameters at 
hand. Even increasing the effective temperature did not 
raise T(0 ++ ) appreciably, because then the ionization pa- 
rameter had to be lowered in order to maintain [O III] 
A5007/[O n] A3727 at its observed level, and this resulted 
in a greater H° abundance, thus enhancing Lya excita- 
tion. Dufour et al. then proposed a composite model, in 
which the [O in] line would be mainly produced around 
high temperature stars (T e // > 38000 K) and the [O II] 
line would be mainly emitted around stars of lower T e f f 
(< 37000 K). Alternatively, one could have, around a star 
of T e ff < 45000 K, a high ionization component emitting 
most of the [O in] and a low ionization component emit- 
ting most of the [O n]. Since then, the HST images (Hunter 
& Thronson 1995, Meurer et al. 1995, Dufour et al. 1996, 
De Mello et al. 1998) have revealed that the NW region 
appears like a shell of ionized gas about 5" in diameter, 
encircling a dense star cluster. Thus the geometries pro- 
posed by Dufour et al. (1988), although quite reasonable 
a priori, do not seem to apply to the case under study. 

Campbell (1990), using Ferland's photoionization code 
CLOUDY, constructed models to fit the spectral obser- 
vations of Lequeux et al. (1979) obtained through a slit 
of 3". 8 x 12". 4. These observations were giving a [O III] 
A4363/5007 ratio of (3.75 ± 0.35) x 10" 2 . With a con- 
stant density, ionization bounded spherical model and a 
LTE Kurucz stellar atmosphere with metallicity 1/10 so- 
lar, in which the adjustable parameters were O/H, T e ff, n 
and e, Campbell obtained a best fit model that had [O III] 
A4363/5007 = 3.07 10~ 2 , i.e. much lower than the value 
she aimed at reproducing. She then proposed a density 
gradient model, in which the inner regions had a density 
n > 10 5 cm -3 , so as to induce collisional deexcitation 
of [O in] A5007. Applying standard abundance deriva- 
tion techniques to this model yields an oxygen abundance 
higher by 70% than the input value. This led Campbell to 
conclude that I Zwl8 was not as oxygen poor as previously 
thought. The density gradient model of Campbell (1990) 
can be checked directly using the density sensitive [Ar iv] 
A4741/4713 ratio. The only observations giving this line 
ratio are those of Legrand et al. (1997), and they indi- 
cate a density in the Ar +++ region lower than 100 cm~ 3 . 
Direct images with the HST do not support Campbell's 
density gradient model either, since, as stated above, the 
appearance of the H n region is that of a shell surrounding 
the excitation stars. 

Stevenson et al. (1993), using a more recent version 
of CLOUDY, constructed a spherical, ionization bounded 



1 A summary of various measurements of [O ill] A4363/5007 
and corresponding electron temperatures is shown in Fig. 2. 



constant density photoionization model to fit their own 
data. They used as an input an extrapolation of the Ku- 
rucz LTE model atmospheres. Their modelling procedure 
was very similar to that of Campbell (1990) for her con- 
stant density model. Their best fit model had O/H = 
1.90 10~ 5 and returned [O ill] A4363/5007 = 2.79 10~ 2 , 
to be compared to their observed value of (3.21 ± 0.42) x 
10~ 2 . 

What complicates the discussion of the three studies 
above is that they use different codes with probably dif- 
ferent atomic data, and they aim at fitting different sets of 
observations. Nevertheless, it is clear that all those models 
have difficulties in reproducing the high [O in] A4363/5007 
observed. They have other weak points, as noted by their 
authors. For example, Dufour et al. (1988) and Stevenson 
et al. (1993) comment on the unsatisfactory fitting of the 
sulfur lines. However, the atomic data concerning sulfur 
are far less well established than those concerning oxygen, 
therefore the discrepancies are not necessarily meaningful. 
Besides, it is not surprising that, with a simple density 
structure, one does not reproduce perfectly at the same 
time the [O m]/[0 n] and [S in] /[S n] ratios. 

The most important defect shared by the three models 
just discussed is that they predict no He n A4686 emission. 
This is simply due to the fact that they used an inadequate 
input stellar radiation field. With the presently available 
stellar population synthesis models for the exciting stars 
of giant H n regions which make use of more realistic 
model atmospheres (Schaerer & Vacca 1998), and espe- 
cially models that are relevant for the Wolf-Rayet stages 
of massive stars, it is interesting to reanalyze the prob- 
lem. Using simple photon counting arguments, De Mello 
at al. (1998) have already shown that a starburst with a 
Salpeter initial mass function and an upper mass limit of 
150 M Q could reproduce the equivalent width of the Wolf- 
Rayet features and of the narrow He II A4686 emission line 
in I Zw 18. 

It is therefore interesting, using the emergent radia- 
tion field from such a synthetic stellar population, to see 
whether one can better reproduce the [O in] A4363/5007 
ratio observed in I Zw 18, with a model that is more com- 
patible with the density structure constrained by the HST 
images. 



3. Our model fitting methodology 

3.1. Computational tools and input parameters 

As in the previous studies, we concentrate on the so-called 
NW component, seen in the top of Fig. 1, which shows the 
WFPC2 Ha image of the HST (cf. Fig. 1 of De Mello et 
al. 1998). Throughout the paper, we adopt a distance to 
I Zw 18 of 10 Mpc, assuming H Q = 75 km s _1 Mpc -1 , as 
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in many studies (Hunter & Thronson 1995, Martin 1996, 
van Zee et al. 1998) |. 

3.1.1. The stellar population 

We use the same model for the stellar population as de- 
scribed in De Mello et al. (1998). It is provided by a evolu- 
tionary population synthesis code using stellar tracks com- 
puted with the Geneva code at the appropriate metallicity 
(1/50 Z Q ). The stellar atmospheres used are spherically 
expanding non-LTE models for WR stars (Schmutz et al. 
1992) and O stars {CoStar models at Z = 0.004, Schaerer 
& de Koter 1997), and Kurucz models at [Fe/H]=-1.5 for 
the remainder. More details can be found in De Mello 
et al. (1998) and Schaerer & Vacca (1998). We assume 
an instantaneous burst of star formation, with an up- 
per mass limit of 150 M and a lower mass limit of 0.8 
M . Since all observational quantities considered here de- 
pend only on the properties of massive stars, the choice 
for Miow has no influence for the results of this paper. 
It merely serves as an absolute normalisation. The total 
initial mass of the stars is adjusted in such a way that, 
at a distance of 10 Mpc, the flux at 3327 A is equal to 
1.7 10~ 15 erg s -1 cm~ 2 A -1 , the value measured in the 
flux calibrated WFPC2 F336W image of De Mello et al. 
(1998) within a circle of 2.5" radius centered on the NW 
region (see Fig. 1). This flux is dominated by the latest 
generation of stars in I Zw 18, so that our normalization is 
hardly sensitive to the previous star formation history in 
the NW region of IZwl8. It yields a total stellar mass of 
8.7 10 4 M , at a distance of 10 Mpc. 

Actually, most of the flux comes from a region much 
smaller in size, and our photoionization modelling is made 
with the ionizing cluster located at the center of the nebula 
and assuming that its spatial extension is negligible. 

We consider that the observed ultraviolet flux is only 
negligibly affected by extinction [J For an extinction 
C(H/3) of 0.04, such as estimated by Izotov & Thuan 
(1998), the corrected flux would be only about 10% larger, 
which is insignificant in our problem. Other observers give 
values of C(H/3) ranging between 0. and 0.2. If C(H/3) 
were as large as 0.20, as estimated by Vflchez & Iglesias- 
Paramo (1998) and some other observers, the true stellar 
flux would be a factor two higher. However, all the deter- 
minations of C(H/3), except the one by Izotov & Thuan 
(1998), do not take into account the underlying stellar ab- 

2 Izotov et al. (1999) have submitted a paper suggesting a 
distance of 20 Mpc to I Zw 18. Should this be the case, the con- 
clusions of our paper that are linked to the ionization structure 
and the temperature of the nebula would hardly be changed. 
The total mass of the ionized gas would be larger, roughly by 
a factor 2 3 . 

3 A direct fitting of the ultraviolet stellar continuum by pop- 
ulation synthesis models, of which we became aware after the 
paper had been submitted, yields C(H/3) < 0.06 (Mas-Hesse & 
Kunth 1999). 



sorption at H/3 and therefore overestimate the reddening. 
A further cause of overestimation of C(H/3), which applies 
also to the work of Izotov & Thuan (1998), is that the 
intrinsic Ha/H/3 ratio is assumed to be the case B recom- 
bination value, while collisional excitation of Ha is not 
negligible in the case of I Zw 18 as noted by Davidson & 
Kinman (1985). We will come back to this below. 

3.1.2. The nebula 

The photoionization computations are performed with the 
code PHOTO using the atomic data listed in Stasihska 
& Leitherer (1996). The code assumes spherical geome- 
try, with a central ionizing source. The diffuse radiation 
is treated assuming that all the photons are emitted out- 
wards in a solid angle of 27r, and the transfer of the res- 
onant photons of hydrogen and helium is computed with 
the same outward only approximation, but multiplying 
the photo-absorption cross-section by an appropriate fac- 
tor to account for the increased path length due to scat- 
tering (Adams 1975). 

The nebular abundances used in the computations are 
those we derived from the spectra of Izotov and Thuan 
(1998) for the NW component of IZwl8, with the same 
atomic data as used in the photoionization code. For he- 
lium, however, we adopted the abundance derived by Izo- 
tov & Thuan (1998) for the SE component, as stellar 
absorption contaminates the neutral helium lines in the 
NW component. The nominal value of the temperature 
derived from [O in] A4363/5007 is 19800 K. This value 
was used to compute the ionic abundances of all the ions 
except + , N + and S + , for which a value of 15000 K was 
adopted (this is the typical value returned by our pho- 
toionization models for IZwl8). The electron density de- 
duced from [S n] A6731 /6717 is 140 cm" 3 , and this density 
was adopted in the computation of the ionic abundances 
of all species. The ionization correction factors to compute 
the total element abundances were those of Kingsburgh & 
Barlow (1994), which are based on photoionization mod- 
els of planetary nebulae and are also suitable for H II re- 
gions. They give slightly smaller oxygen abundances (by 
a few %) than the traditional ionization correction factors 
which assume that the +++ region is coextensive with 
the He ++ (we did not iterate on the ionization correc- 
tion factors after our photoionization model analysis since 
this would have not changed any of the conclusions drawn 
in this paper). The carbon abundance used in the com- 
putations follows from the C/O ratio derived by Garnett 
et al. (1997) from HST observations of IZwl8. The abun- 
dances of the elements not constrained by the observations 
(Mg, Si) and (CI, Fe) have been fixed to 10~ 7 and 10~ 8 
respectively. Table 1 presents the abundance set used in 
all the computations presented in the paper. As already 
noted by previous authors, at the metallicity of I Zw 18, 
the heavy elements (i.e. all the elements except hydrogen 
and helium) play a secondary role in the thermal balance. 
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Table 1. Input abundances relative to hydrogen for the 
models of I Zw 18 (by number). 



He 


7.60 10" 


-'2 


C 


3.03 10" 


■6 


N 


3.89 10" 


-7 





1.32 10" 


■5 


Ne 


2.28 10" 


■6 


S 


3.72 10" 


■7 


Ar 


9.13 10" 


■8 



Their role in the absorption of ionizing photons is com- 
pletely negligible. Any change of abundances, even that 
of helium, compatible with the observed intensities of the 
strong lines, will result in a very small change in the elec- 
tron temperature, and we have checked that the effect they 
will induce in the predicted spectra are small compared to 
the effects discussed below. 

We do not include dust in the computations. While 
it is known that, in general, dust mixed with the ionized 
gas may absorb some of the ionizing photons, and con- 
tribute to the energy balance of the gas by photoelectric 
heating and collisional cooling (e.g. Baldwin et al. 1991, 
Borkowski & Harrington 1991, Shields & Kennicutt 1995), 
the expected effect in I Zw 18 is negligible, since the dust- 
to-gas ratio is believed to be small at such metallicities 
(cf. Lisenfeld & Ferrara 1998). 

The case of IZwl8 is thus very interesting for pho- 
toionization modelling, since due to the very low metal- 
licity of this object, the number of unconstrained relevant 
parameters is minimal. 

3.2. Fitting the observational constraints 

In judging the value of our photoionization models, we 
do not follow the common procedure of producing a table 
of intensities relative to H/3 to be compared to the ob- 
servations. A good photoionization model is not only one 
which reproduces the observed line ratios within the un- 
certainties. It must also satisfy other criteria, like being 
compatible with what is known from the distribution of 
the ionized gas, and what is known of the ionizing stars 
themselves. On the other hand, many line ratios are not 
at all indicative of the quality of a photoionization model. 
For example, obviously, two lines arising from the same 
atomic level like [O in] A5007 and [O in] A4959 have in- 
tensity ratios that depend only on the respective transition 
probabilities. In H n regions, the ratio of hydrogen Balmer 
lines (if case B applies) is little dependent on the physical 
conditions in the ionized gas, and this is why it can be 
used to determine the reddening. The ratios of the inten- 
sities of neutral helium lines do depend somewhat on the 
electron density distribution and on selective absorption 
by dust of pseudo-resonant photons, (Clegg & Harrington 
1989, Kingdon & Ferland 1995), and these are introduced 



in photoionization models. In the case of the NW com- 
ponent of I Zw 18, the observed neutral helium lines are 
affected by absorption from stars or interstellar sodium 
(Izotov & Thuan 1998, Vflchez & Iglesias-Paramo 1998), 
and cannot be easily used as constraints for photoioniza- 
tion models. 

Generally speaking, once line ratios indicative of 
the electron temperature (like [O in] A4363/5007, 
[N n] A5755/6584), of the electron density (like [S n] 
A6731/6717, [Ar iv] A4741/4713) and of the global ion- 
ization structure (like [O ill] A5007/[O n] A3727 or [S in] 
A9532/[S Ii] A6725) have been fitted, the ratios of all the 
strong lines with respect to H/3 are necessarily reproduced 
by a photoionization model whose input abundances were 
obtained from the observations. The only condition is that 
the atomic data to derive the abundances and to compute 
the models should be the same. Problems may arise only 
if the empirical ionization factors are different from the 
ones given by the model, or if there is insufficient infor- 
mation on the distribution of the electron temperature or 
density inside the nebula (in the case of I Zw 18 no direct 
information is available on the temperature in the low ion- 
ization zone, but we adopted a value inspired by the mod- 
els). Therefore, intensity ratios such as [O III] A5007/H/3, 
[Ne in] A3869/H/3, [N n] A6584/H/3, [Ar in] A7135/H/3 or 
C in] A1909/H/3 are not a measure of the quality of the 
photoionization model. 

To judge whether a photoionization model is accept- 
able, one must work with outputs that are significantly 
affected by the physical processes on which the photoion- 
ization model is based, i.e. the transfer of the ionizing 
radiation, the processes determining the ionization equi- 
librium of the various atomic species and the thermal bal- 
ance of the gas. Table 2 lists the quantities that can be 
used in the case of IZwl8, given the observational in- 
formation we have on the object. The value of the H/3 
flux is derived from the Ha flux measured in a circle of 
radius 9—2.5" (shown in Fig. 1), assuming C(H/3) = 0. 
The line ratios He n A4686/H/3, [O in] A4363/5007, [S n] 
A6731/6717, [O in] A5007/[O n] A3727, [S m] A6312/[S n] 
A6725 and [O i] A6300/H/3 are the values observed by Izo- 
tov & Thuan (1998) in the rectangular aperture whose 
approximate position is shown in Fig. 1. It is important 
to define in advance the tolerance we accept for the differ- 
ence between our model predictions and the observations. 
This must take into account both the uncertainty in the 
observational data, the fact that the spectra were taken 
through an aperture not encompassing the whole nebula, 
and the fact that the nebula does not have a perfect, spher- 
ical symmetric structure. This latter aspect is, of course, 
difficult to quantify, and the numbers given in Column 3 of 
Table 2 are to be regarded rather as guidelines. In Column 
4, we indicate which is the dominant factor determining 
the adopted tolerance : the signal-to-noise, or the geom- 
etry. For example, such ratios as He n/H/3, [O ni]/[0 Ii], 
[S in] /[S ii] or [O i]/H/3 are obviously more dependent 
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on geometrical effects than [O in] A4363/5007. Note that, 
even for that ratio, the tolerance given in Table 2 is larger 
than the uncertainty quoted by Izotov & Thuan (1998). 
The reason is that the many observations of the NW com- 
ponent of I Zw 18 made over the years, with different tele- 
scopes, detectors, and apertures, yield distinct values for 
this ratio, as shown in Fig. 2. In view of this figure, a tol- 
erance of 10% with respect with the [O ill] A4363/5007 
ratio measured by Izotov & Thuan (1998) seems reason- 
able. The status of the [S n] A6731/6717 ratio is somewhat 
different. It indicates the average electron density in the 
zone emitting [S n] . This is very close to an input param- 
eter, since photoionization models are built with a given 
density structure. However, because the density deduced 
from [S n] A6731/6717 is not the total hydrogen density 
but the electron density in the region emitting [S n] , and 
because the density is not necessarily uniform, it is impor- 
tant to check that the model returns an [S n] A6731/6717 
value that is compatible with the observations. For the 
total H/3 flux, we accept models giving F(H/3) larger than 
the observed value, on account of the fact that the cov- 
erage factor of the ionizing source by the nebula may be 
smaller than one. 

Thus, in the following, we compute photoionization 
models with the input parameters as defined in Section 
3.1, and see how they compare with the constraints speci- 
fied in Table 2. We will not examine the effects of varying 
the elemental abundances, since, as mentioned above, they 
are negligible in our problem. Uncertainties in the mea- 
sured stellar flux have only a small impact on our models, 
and are therefore not discussed here. Similarly, we discard 
the effects of an error in the distance d to I Zw 18. These 
are not crucial on the ionization structure of a model de- 
signed to fit the observed flux at 3327A, since the mean 
ionization parameter varies roughly like <i 2 / 3 . What we 
mainly want to see is whether, with our present knowl- 
edge, we can satisfactorily explain the observed properties 
of I Zw 18. As will be seen, the gas density distribution 
plays an important role. 

4. Climbing the ladder of sophistication 

4--1- The ionizing radiation field 

Before turning to proper photoionization modelling, it is 
worthwhile examining the gross properties of the ionizing 
radiation field of the synthetic stellar population model 
we are using, and compare it to single star model atmo- 
spheres. Two quantities are particularly relevant. One is 
<3(Hc + )/<3(H ), the ratio of the number of photons above 
54.4 and 13.6 eV emitted by the ionizing source. This ratio 
allows one to estimate the He n A4686/H/3 that would be 
observed in a surrounding nebula, by using simple conser- 
vation arguments leading to the formula: He II A4686/H/3 
= 2.14 Q(He + )/Q(H°) (taking the case B recombination 
coefficients given in Osterbrock 1989). As is known, this 
expression is valid if the nebula is ionization bounded and 



the observations pertain to the whole volume. It is less 
commonly realized that it also assumes the average tem- 
perature in the Hc ++ region to be the same as in the 
entire H n region. This is may be far from true, as will 
be shown below, so a correction should account for that. 
Another assumption is that all the photons above 54.4 eV 
are absorbed by He + ions. This is not what happens in 
objects with a low ionization parameter. There, the resid- 
ual neutral hydrogen particles are sufficiently numerous to 
compete with He + . In such a case, the expression above 
gives an upper limit to the nebular He II A4686/H/3 ratio. 
In spite of these difficulties, Q(Rc + )/Q(R°) remains a use- 
ful quantity to estimate the intensity of the nebular He n 
A4686 line. Fig. 3a shows the variation of Q(Hc+)/Q(H°) 
as a function of starburst age for the synthetic model pop- 
ulation we arc considering. As already stated in De Mello 
et al. (1998), the strength of the He II nebular line in 
IZw 18 indicates a starburst between 2.9 and 3.2 Myr. 

Another important ratio is Q(He°)/Q(H°), sometimes 
referred to as the hardness parameter of the ionizing radi- 
ation field. It provides a qualitative measure of the heat- 
ing power of the stars. We have represented this quantity 
in Fig. 3c. We see that, as the starburst ages, its heat- 
ing power gradually declines, and shows only a very mild 
bump at ages around 3 Myr, where the Wolf-Rayet stars 
are present. As we will show below this modest increase 
of the heating power is not sufficient to explain the high 
electron temperature observed in IZw 18. 

For comparison, we show in Figs. 3b and d respectively, 
the values of Q(He+)/Q(H°) and Q(He°)/Q(H°) as a func- 
tion of the stellar effective temperature for the LTE model 
atmospheres of Kurucz (1991), the CoStar model atmo- 
spheres corresponding to main sequence stars (Schaerer & 
dc Koter 1997) and for the model atmospheres for Wolf- 
Rayet stars of Schmutz et al. (1992). The CoStar mod- 
els show an increased He + ionizing flux compared to Ku- 
rucz models which have a negligible flux even for very low 
metallicities ([Fe/H]= —1.5). The reasons for this differ- 
ence have been discussed in Schaerer & de Koter (1997). In 
addition to the T e g dependence, Q(He + ) from WR models 
depend strongly on the wind density. He + ionizing photons 
are only predicted by models with sufficiently thin winds 
(cf. Schmutz et al. 1992). Figure 3d shows the increase of 
the hardness of the radiation field, at a given T e s, between 
the spherically expanding non-LTE models for O and WR 
stars and the traditional Kurucz models (see discussion in 
Schaerer & de Koter 1997). This provides a greater heat- 
ing power which, as will be shown later, is however still 
insufficient to explain the observations. 

4-2. IZw 18 as a uniform sphere 

We start with the easiest and most commonly used ge- 
ometry in photoionization modelling: a sphere uniformly 
filled with gas at constant density, occupying a fraction e 
of the whole nebular volume. The free parameters of the 
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Table 2. Observables that a photoionization model of IZwl8 should fit (references in Sect. 3.2). 



Quantity 


Value 


Tolerance 


Major source 


of uncertainty 


Symbol in Figs. 4-6 


F(H/3) [erg cm" 2 s" 1 ] 


4.0 10 -14 


+ 0.5 dex 


geometry (see 


text) 


circle 


angular radius 6 [arc sec ] 


2.5 


± 0.08 dex 


geometry 




cross 


He ii A4686/H/3 


0.034 


± 0.2 dex 


geometry (see 


text) 


square 


[0 in] A4363/5007 


3.28 10~ 2 


± 0.04 dex 


S/N 




open triangle 


[S n] A6731/6717 


1.3 


± 0.04 dex 


S/N 




diamond 


[0 in] A5007/[O n] A3727 


6.82 


± 0.1 dex 


geometry 




filled triangle 


[0 i] A6300/H/3 


0.007 


± 0.3 dex 


geometry 




plus 


[S m] A6312/[S n] A6725 


0.173 


± 0.2 dex 


geometry 




asterisk 



models are then only the age of the starburst, the gas den- 
sity and the filling factor. Each model is computed start- 
ing from the center, and the computations are stopped 
either when the [O m]/[0 n] ratio has reached the ob- 
served value given in Table 2, or when the gas becomes 
neutral. In other words, we examine also models that are 
not ionization bounded, in contrast to previous studies. 

Figure 4 shows our diagnostic diagram for a series of 
models having a density n= 100 cm -3 and a filling fac- 
tor e=0.01. The left panel shows the computed values of 
log F(H/3) + 15 (open circles), log He n A4686/H/3 + 2 
(squares), angular radius 9 (crosses), [O in] A4363/5007 
x 100 (open triangles), [S ii] A6731/6717 (diamonds), log 
([O m]/[0 ii]) (black triangles), log ([S m]/[S ii]) (as- 
terisks) and log [O i]/H/3 +3 (plus) as a function of the 
starburst age. The black circles correspond to the value 
of log F(H/3) + 15 that the nebula would have if it were 
ionization bounded. Thus, by comparing the positions of 
an open circle and a black circle, at a given abscissa, one 
can immediately see whether the model is density bounded 
and how much diffuse Ha or H/3 emission is expected to be 
emitted outside the main body of the nebula. In the right 
panel, the observed values are represented on the same 
vertical scale and with the same symbols as the model 
predictions. The tolerances listed in Table 2 are repre- 
sented as vertical error bars (the horizontal displacement 
of the symbols has no particular meaning). We readily see 
that the age of the starburst is important only for the He II 
A4686 line, the other quantities varying very little for ages 
2.7-3.4 Myr. Therefore, for the following runs of models, 
we adopt an age of 3.1 Myr. In principle, one can always 
adjust the age for the model to reproduce the observed 
He II A4686/H/3 ratio exactly. 

Figure 5 shows the same sort of diagnostic diagram as 
Fig. 4 for a series of models with n = 100 cm~ 3 and vary- 
ing filling factor. For a filling factor around 0.1 or larger, 
with the adopted electron density, the model is ionization 
bounded, and its [O m]/[0 n] is larger than observed. For 
filling factors smaller than that, the gas distribution is 
more extended, so that the general ionization level drops. 
The observed [O m]/[0 ii] can then only be reproduced 
for a density bounded model. In such a case, the H/3 radia- 



tion actually produced by the nebula is smaller than if all 
the ionizing photons were absorbed in the nebula. A filling 
factor of 0.002 - 0.05 gives values of [O m]/[0 n], F(H/3) 
and 9 in agreement with the observations. But such mod- 
els give [O in] A4363/5007 too small compared with the 
observations, and [O i]/H/3 below the observed value by 
nearly two orders of magnitude. It is interesting, though, 
to understand the qualitative behavior of these line ra- 
tios as e decreases. [O i]/H/3 decreases because the model 
becomes more and more density bounded in order to re- 
produce the observed [O m]/[0 n] and levels off at e = 
0.1, because the ionization parameter of the model is then 
so small that the [O i] is gradually emitted by residual 
neutral oxygen in the main body of the nebula and not in 
the outskirts. [O in] A4363/5007 decreases as e decreases, 
because of the increasing proportion of La cooling as the 
ionization parameter drops. 

One can build other series of models with different val- 
ues of n that arc still compatible with the observed [S n] 
A6731/6717. Qualitatively, their behavior is the same and 
no acceptable solution is found. 

Interestingly, Fig. 5 shows that models with e > 0.1 
have [O in] A4363/5007 marginally compatible with the 
observations ( [O in] A4363/5007 = 3.03 10~ 2 for e = 0.1), 
but such models have too large [O m]/[0 n](> 10 com- 
pared of the observed value 6.8) and too small angular ra- 
dius (< 1.6" instead of the observed value 2.5"). Note, by 
the way, that such models, being optically thick, return a 
rather large [O i]/H/3, actually close to the observed value, 
and a [S m]/[S ii] compatible with the observations. How- 
ever, we do not use [S m]/[S ii] as a primary criterion to 
judge the validity of a model, since experience with pho- 
toionization modelling of planetary nebulae shows that it 
is difficult to reproduce at the same time the sulfur and 
the oxygen ionization structure of a given object, and, in 
principle, one expects the atomic data for oxygen to be 
more reliable than those for sulfur. The strongest argu- 
ment against models with e > 0.1 is their angular size, 
which is definitely too small compared with the observa- 
tions. This remains true even when considering a reason- 
able error on the distance since, with the condition that we 
impose on the flux at 3327 A to be preserved, the angular 
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radius of a model goes roughly like d 1 / 3 . This illustrates 
the importance of taking into account other parameters in 
addition to line ratios to accept or reject a photoionization 
model. 

4-3. I Zw 18 as a spherical shell 

The series of models presented above had mainly a ped- 
agogical value, but they are obviously incompatible with 
the observed morphology in Ha. The next step is to con- 
sider a model consisting of a hollow, spherical shell of con- 
stant density, similar to the one constructed by Garcia- 
Vargas et al. (1997) for NGC 7714 for example. In such a 
case, there is an additional free parameter, Ri ni the ra- 
dius of the inner boundary of the shell. It is fixed, more or 
less, by the appearance of the Ha image. Figure 6 shows a 
diagnostic diagram for a series of models with R in — 2.25 
10 20 cm (corresponding to an angular radius of 1.5"), n = 
100 cm~ 3 and varying e. The qualitative behavior is simi- 
lar to that seen for the uniform sphere models presented in 
Fig. 1, but the [O in] A4363/5007 ratio is now even lower 
(it never exceeds 2.5 10~ 2 in this series). This is because 
of the enhanced role of La cooling, which is strong in all 
the parts of the nebula, while for the full sphere model, in 
the zone close to the star, the ionization parameter is very 
high and consequently the population of neutral hydrogen 
very small. 

Apart from the [O in] A4363/5007 problem, models 
with e = 0.002 - 0.05 are satisfactory as concerns the 
main diagnostics ([O m]/[0 n], F(H/3) and 6). The mod- 
els become progressively density bounded towards smaller 
values of e < 0.02, meaning that there is a leakage of ioniz- 
ing photons. From Fig. 6, one sees that these photons are 
enough to produce an Ha emission in an extended diffuse 
envelope that is at least comparable in strength to the 
total emission from the dense shell. This is in agreement 
with Dufour & Hester's (1990) ground-based observation 
of extended Ha emission surrounding the main body of 
star formation. 

4-4- Other geometries 

Closer inspection of the HST Ha image shows that the 
gas emission is incompatible with a spherical bubble. This 
is illustrated in Fig. 7, where the observed cumulative 
surface brightness profile within radius r (dashed line) 
is compared to the expected profiles for constant density 
spherical shells of various inner radii (solid lines). The 
theoretical profiles are obtained assuming that the tem- 
perature is uniform in the gas, but taking into account 
a reasonable temperature gradient in the model hardly 
changes the picture. Clearly, the observed profile is not 
compatible with a spherical density distribution. The col- 
umn density of emitting matter in the central zone of the 
image is too small. One must either have an incomplete 
shell with some matter stripped off from the poles, or even 



a more extreme morphology like an diffuse axisymmetric 
body with a dense ringlike equator seen face on. Such ge- 
ometries are actually common among planetary nebulae 
(Corradi & Schwartz 1995) and nebulae surrounding lu- 
minous blue variables (Nota et al. 1995), being probably 
the result of the interaction of an aspherical stellar wind 
from the central stars (Mcllcma 1995, Frank et al. 1998) 
and are also suggested to exist in superbubbles and super- 
giant shells and to give rise to the blow-out phenomenon 
(Mac Low et al. 1989, Wang & Helfand 1991, Tcnorio- 
Tagle et al. 1997, Oey & Smedley 1998, Martin 1998). 

Does the consideration of such a geometry help in solv- 
ing the [O in] A4363/5007 problem? In the case of a spher- 
ical shell with some matter lacking at the poles, our com- 
putations overestimate the role of the diffuse ionizing ra- 
diation, which is supposed to come from a complete shell. 
Since the heating power of the diffuse ionizing radiation is 
smaller than that of the stellar radiation, one may be un- 
derestimating the electron temperature. As a test, we have 
run a model switching off completely the diffuse radiation, 
so as to maximize the electron temperature. The effect was 
to increase the ++ temperature by only 200 K. In the 
case of a ring with some diffuse matter seen in projection 
inside the ring, the gas lying close to the stars would be 
at a higher electron temperature than the matter of the 
ring, and one expects that an adequate combination of the 
parameters describing the gas distribution interior to and 
inside of the ring might reproduce the [O ill] A4363/5007 
measured in the aperture shown in Fig. 1. However, we 
notice that the region of high [O in] A4363/5007 is much 
larger than that. It extends over almost 20" ( Vflchez 
& Iglesias-Paramo 1998) and [O in] A4363/5007 is nei- 
ther correlated with [O m]/[0 il] nor with the Ha surface 
brightness. Therefore, one cannot explain the high [O in] 
A4363 by the emitting matter being close to the ionizing 
stars. In passing we note that the observations of Vflchez 
& Iglesias-Paramo (1998) show the nebular He n A4686 
emission to be extended as well (although not as much as 
[O in] A4363). If this emission is due to photoionization 
by the central star cluster, as modelled in this paper, this 
means that the Ha ring is porous, since He n A4686 emis- 
sion necessarily comes from a region separated from the 
stars by only a small amount of intervening matter. 

In summary, we must conclude that, if we take into 
account all the available information on the structure of 
I Zw 18, we are not able to explain the high observed [O in] 
A4363/5007 ratio with our photoionization models. In our 
best models, [O in] A4363/5007 is below the nominal value 
of Izotov & Thuan (1998) by 25 - 35%. 

4-5. Back to the model assumptions 

Can our lack of success be attributed to an improper de- 
scription of the stellar radiation field? After all, we know 
little about the validity of stellar model atmospheres in 
the Lyman continuum (see discussion in Schaerer 1998). 
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Direct measurements of the EUV flux of early B stars re- 
vealed an important EUV excess (up to <~ 1.5 dex) with re- 
spect to plane parallel model atmospheres (Cassinelli et al. 
1995, 1996), whose origin has been discussed by Najarro et 
al. (1996), Schaerer & de Koter (1997) and Aufdenberg et 
al. (1998). For O stars a similar excess in the total Lyman 
continuum output is, however, excluded from considera- 
tions of their bolomctric luminosity and measurements of 
H II region luminosities (Oey & Kennicutt 1997, Schaerer 
1998). The hardness of the radiation field, which is cru- 
cial for the heating of the nebula, is more difficult to test. 
Some constraints can be obtained by comparing the line 
emission of nebulae surrounding hot stars with the results 
of photoionization models (Esteban et al. 1993, Peha et 
al. 1998, Crowther et al. 1999), but this is a difficult task, 
considering that the nebular emission depends also on its 
geometry. Although the hardness predicted by the CoStar 

stars models permits to build grids of photoionization 
models that seem to explain the observations of Galac- 
tic and LMC H n regions (Stasihska & Schaerer 1997), 
the constraints are not sufficient to prove or disprove the 
models. 

To check the effect of a harder radiation field, we have 
run a series of models where the radiation field above 
24.6 eV was arbitrarily multiplied by a factor 3 (raising 
the value of Q(Hc°)/Q(H°) from 0.33 to 0.59, correspond- 
ing to T cff from ^ 40000 K to - 100 000 K or a blackbody 
of the same temperature (cf. Fig. 3d). This drastic hard- 
ening of the radiation field resulted in an increase of [O in] 
A4363/5007 of only 10%. It is only by assuming a black- 
body radiation of 300 000 K (which has Q(He°)/Q(H°) = 
0.9) that one approaches the observed [O in] A4363/5007. 
A model similar to those presented in Fig. 6 but with such 
a radiation field gives [O in] A4363/5007 =3.15 10~ 2 . But 
is has a He II A4686/H/3 ratio of 0.53, which is completely 
ruled out by the observations. Of course, a blackbody is 
probably not the best representation for the spectral en- 
ergy distribution of the radiation emitted by a very hot 
body, but in order to explain the emission line spectrum of 

1 Zw 18 by stars, one has to assume an ionizing radiation 
strongly enhanced at energies between 20 - 54 eV, but 
not above 54.4 eV, compared to the model of De Mello et 
al. (1998). If this is realistic cannot be said at the present 
time. 

We have also checked the effect of heating by addi- 
tional X-rays that would be emitted due to the interac- 
tion of the stellar winds with ambient matter (see Mar- 
tin 1996 for such X-ray observations), by simply adding a 
bremsstrahlung spectrum at T = 10 6 or T = 10 7 K to the 
radiation from the starburst model. As expected, the ef- 
fect on [O ill] A4363/5007 was negligible, since the X-rays 
are mostly absorbed in the low ionization regions (they do 
raise the temperature in the + zone to T e ~ 16 000 K). 

As already commented by previous authors, changing 
the elemental abundances does not improve the situation. 
Actually, even by putting all the elements heavier than he- 



lium to nearly zero abundance, [O in] A4363/5007 is raised 
by only 7%. Varying the helium abundance in reasonable 
limits does not change the problem. 

The neglect of dust is not expected to be responsible 
for the too low [O III] A4363/5007 we find in our models 
for I Zw 18. Gas heating by photoemission from grains 
can contribute by as much as 20% to the electron thermal 
balance when the dust-to-gas ratio is similar to that found 
in the Orion nebula. But, as discussed by Baldwin et al. 
(1991), it is effective close to the ionizing source where 
dust provides most of the opacity. Besides, the proportion 
of dust in I Zw 18 is expected to be small, given the low 
metallicity. Extrapolating the relation found by Liscnfeld 
& Ferrara (1998) between the dust-to-gas ratio and the 
oxygen abundances in dwarf galaxies to the metallicity of 
I Zw 18 yields a dust-to-gass mass ratio 2 to 3 orders of 
magnitudes smaller than in the solar vicinity. 

It remains to examine the ingredients of the photoion- 
ization code. We first stress that the [O m] A4363/5007 
problem in I Zw 18 has been encountered by various au- 
thors using different codes, even if the observational ma- 
terial and the points of emphasis in the discussion were 
not the same in all the studies. As a further test, we com- 
pared the same model for I Zw 18 run by CLOUDY 90 
and by PHOTO and the difference in the predicted value 
of [O in] A4363/5007 was only 5%. One must therefore in- 
criminate something that is wrongly treated in the same 
way in many codes. One possibility that comes to mind 
is the diffuse ionizing radiation, which is treated by some 
kind of outward only approximation in all codes used to 
model I Zw 18. However, we do not expect that an ac- 
curate treatment of the diffuse ionizing radiation would 
solve the problem. Indeed, comparison of codes that treat 
more accurately the ionizing radiation with those using an 
outward only approximation shows only a negligible dif- 
ference in the electron temperature (Ferland et al. 1996). 
Besides, as we have shown, even quenching the diffuse ion- 
izing radiation does not solve the problem. 

Finally, one can also question the atomic data. The 
most relevant ones here are those governing the emissiv- 
ities of the observed [O III] lines and the H I collision 
strengths. The magnitude of the discrepancy we wish to 
solve would require modifications of the collision strengths 
or transition probabilities for [O III] of about 25%. This 
is much larger than the expected uncertainties and the 
differences between the results of different computations 
for this ion (see discussion in Lennon and Burke 1994 and 
Galavis et al. 1997). Concerning La excitation, dividing 
the collision strength by a factor 2 (which is far above 
any conceivable uncertainty, see e.g. Aggarwal et al. 1991) 
modifies [O in] A4363/5007 only by 2% because La acts 
like a thermostat. 

We are therefore left with the conclusion that the 
[O in] A4363/5007 ratio cannot be explained in the frame- 
work of photoionization models alone. 
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4-6. Condensations and filaments 

Another failure of our photoionization models is that they 
predict too low [O i]/H/3 compared to observations. This 
is a general feature of photoionization models and it is 
often taken as one of the reasons to invoke the presence 
of shocks. However, it is well known that the presence 
of small condensations of filaments enhances the [O i]/H/3 
ratio, by reducing the local ionization parameter. Another 
possibility is to have an intervening filament located at a 
large distance from the ionizing source, whose projected 
surface on the aperture of the observations would be small. 

In order to see under what conditions such models can 
quantitatively account for the observed [O i]/H/3 ratio, we 
have mimicked such a situation by computing a series of 
ionization bounded photoionization models for filaments 
of different densities located at various distances from the 
exciting stars. For simplicity, we assumed that there is no 
intervening matter between the filaments and the ionizing 
source. The models were actually computed for complete 
shells. The radiation coming from a filament can be simply 
obtained by multiplying the flux computed in the model 
by the covering factor / by which the filament is covering 
the source. 

Table 3 presents the ratios with respect to H/3 of the 
[O i] A6300, [O n] A3727, [O in] A5007, [S n] A6717 and 
[S n] A6731 for these models. The first column of the ta- 
ble corresponds to a photoionization model for the main 
body of the nebula with R in = 2.25 10 20 cm (correspond- 
ing to #in = 1-5"), "=100 cm~ 3 , and e=0.01, and density 
bounded so as to obtain the observed [O m]/[0 Ii] (this is 
one of the models of Fig. 6). It can be readily seen that, 
for an intervening filament of density 10 2 cm -3 located at 
a distance of 500 pc from the star cluster, or for condensa- 
tions of density 10 6 cm~ 3 , one can reproduce the observed 
[O i]/H/3 line ratio without strongly affecting the remain- 
ing lines, not even the density sensitive [S Ii] A6731/6717 
ratio, if one assumes a covering factor / of about 0.1. 

This explanation may appear somewhat speculative. 
However, one must be aware that the morphology of the 
ionized gas in I Zw 18 shows filaments in the main shell 
as well as further out and this has been amply discussed 
in the literature (Hunter & Thronson 1995, Dufour et al. 
1995, Martin 1996). Our point is that, even in the frame- 
work of pure photoionization models, if one accounts for a 
density structure suggested directly by the observations, 
the strength of the [O i] line can be easily understood. 

We note that the condensations or filaments that pro- 
duce [O i] are optically thick, and therefore their neu- 
tral counterpart should be seen in H I. Unfortunately, the 
available H I maps of I Zw 18 (van Zee et al. 1998) do not 
have sufficient spatial resolution to reveal such filaments. 
But these authors show that the peak emission in H I and 
Ha coincide (in their paper, the peak emission in Ha ac- 
tually refers to the whole NW component), and find that 
the entire optical system, including the diffuse emission, is 



embedded in a diffuse, irregular and clumpy neutral cloud. 
In such a situation, it is very likely that some clumps or 
filaments, situated in front of the main nebula, and having 
a small covering factor, produce the observed [O i] line. 

By using photoionization models to explain the ob- 
served [O i]/H/3 ratio, one can deduce the presence of den- 
sity fluctuations, even if those are not seen directly. Such 
density fluctuations could then be, more directly, inferred 
from density sensitive ratios of [Fe n] lines, such as seen 
and analyzed in the Orion nebula by Bautista & Pradhan 
(1995). 

4-7. A few properties of IZwl8 deduced from models 

Although we have not built a completely satisfactory pho- 
toionization model reproducing all the relevant observed 
properties of IZwl8, we are not too far from it. We de- 
scribe below some properties of the best models that may 
be of interest for further empirical studies of this object. 

Tables 4 and 5 present the mean ionic fractions and 
mean ionic temperatures (as defined in Stasinska 1990) 
for the different elements considered, in the case of the 
best fit models with a uniform sphere and a spherical shell 
respectively, with n=100 cm -3 , and e=0.01. 

It can be seen that, while both geometries yield similar 
ionic fractions for the main ionization stages, the relative 
populations of the highly charged trace ions are very dif- 
ferent. In the case of the uniform sphere, the proportion 
of +++ , for example, is twice as large as in the shell 
model. Also, the characteristic temperature of ions with 
high ionization potential are much higher in the case of 
the filled sphere, for reasons commented on earlier. As a 
result, the OIV]1394/H/3 ratio is 9.7 10" 3 in the first case 
and 8.6 10 -4 in the second. This line is too weak to be 
measured in I Zw 18, of course, but it is useful to keep this 
example in mind for the study or more metal rich objects. 

It is interesting to point out that, in the case of the 
uniform sphere, the total flux in the He II A4686 line is 
smaller than in the case of the shell model (3.1 10~ 15 vs. 
3.3 10~ 15 erg cm~ 2 s _1 ) despite the fact that the ionic 
fractions of He ++ are similar. This is because the Hc ++ 
region is at a much higher temperature (25 000 K versus 
18000 K). 

Tables 4 and 5 can be used for estimating the ioniza- 
tion correction factors for IZwl8. Caution should, how- 
ever, be applied regarding especially the ionization struc- 
ture predicted for elements of the third row of Mendeleev's 
table. Experience in photoionization modelling of plan- 
etary nebulae, where the observational constraints are 
larger, shows that the ionization structure of these ele- 
ments is rarely satisfactorily reproduced with simple mod- 
els (Howard et al. 1997, Pcna et al., 1998). 

The total ionized mass in the NW component is rel- 
atively well determined, since we know the total num- 
ber of ionizing photons, the radius of the emitting region 
and have an estimate of the mean ionization parameter 
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Main body 


Filaments of various 


n 


Filaments at various distances 


n [cm 3 ] 


10 2 a 


10 4 


10 5 


10 6 


10 2 


10 2 


Oin["] 




1.5 


1.5 


1.5 


20. 


100. 


F(H/3) f ' 


2.02e-13 


3.33e-13 


3.32e-13 


3.35e-13 


3.41e-13 


3.47e-13 


[0 i] A6300 


1.04e-4 


5.28e-2 


1.57e-l 


3.83e-l 


6.90e-2 


2.81e-l 


[0 n] A3727 


2.15e-l 


2.41e-l 


5.39e-2 


1.19e-2 


4.71e-l 


2.85e-l 


[0 in] A4363 


3.46e-2 


3.63e-5 


1.12e-6 


1.04e-7 


1.20e-5 


1.74e-7 


[0 in] A5007 


1.42 


3.76e-3 


1.41e-4 


1.25e-5 


1.41e-3 


3.08e-5 


[S n] A6717 


5.58e-4 


5.20e-2 


5.23e-2 


5.27e-2 


1.87e-l 


5.70e-l 


[S II] A6731. 


4.19e-3 


8.04e-2 


9.16e-2 


1.03e-l 


1.35e-l 


3.96e-l 



a The model for the main body is density bounded (see text) 
b The total H/3 flux is given in erg s _1 for a covering factor f=l. 



through the observed [O m]/[0 n]. Indeed, at a given 
Q(H°), for a constant density sphere with a filling factor 
e, n 2 e is proportional to the cube of the radius, while ne 2 
is proportional to the cube of the ionization parameter. Of 
course, we have just made the point that the NW compo- 
nent of I Zw 18 is not a sphere. Nevertheless, we have an 
order of magnitude estimate, which is of 3. 10 5 M at d 
= 10 Mpc (this estimate varies like d 3 ). 

Finally, it is important to stress, as already mentioned 
above, that Ha is partly excited by collisions. In all our 
models for the main body of the nebula, Ha/H/3 lies be- 
tween 2.95 and 3, while the case B recombination value is 
2.7. This means that the reddening of IZwl8 is smaller 
than the value obtained using the case B recombination 
value at the temperature derived from [O in] A4363/5007. 
If we take the observations of Izotov & Thuan (1998), who 
also correct for underlying stellar absorption, we obtain 
C(H/3)=0. 



5. Summary and conclusions 

We have built photoionization models for the NW compo- 
nent of I Zw 18 using the radiation field from a starburst 
population synthesis at appropriate metallicity (De Mcllo 
et al. 1998) that is consistent with the Wolf-Rayet sig- 
natures seen in the spectra of I Zw 18. The aim was to 
see whether, with a nebular density structure compatible 
with recent HST images, it was possible to explain the 
high [O in] A4363/5007 ratio seen in this object, com- 
monly interpreted as indicative of electron temperature ~ 
20 000K. 

For our photoionization analysis we have focused on 
properties which are relevant and crucial model predic- 
tions. For the observational constraints we have not only 
used line ratios, but also other properties such as the in- 
tegrated stellar flux at 3327 A and the observed angular 
radius of the emitting region as seen by the HST. Care has 
also been taken to include tolerances on model properties 



which may be affected by deviations from various simple 
geometries. 

We have found that [O in] A4363/5007 cannot be ex- 
plained by pure photoionization models, which yield too 
low an electron temperature. We have considered the ef- 
fects due to departure from spherical symmetry indicated 
by the HST images. Indeed these show that the NW com- 
ponent of I Zw 18 is neither a uniform sphere nor a spheri- 
cal shell, but rather a bipolar structure with a dense equa- 
torial ring seen pole on. 

We have discussed the consequences that an inaccu- 
rate description of the stellar ionizing radiation field might 
have on [O in] A4363/5007 as well as additional photoion- 
ization by X-rays. Finally, we have considered possible er- 
rors in the atomic data. 

All these trials were far from solving the electron tem- 
perature problem, raising the [O in] A4363/5007 ratio by 
only a few percent while the discrepancy with the observa- 
tions is on the 30% level. Such a discrepancy means that 
we are missing a heating source whose power may be of 
the same magnitude as that of the stellar ionizing photons. 
It is also possible that the unknown energy source is not 
so powerful, but acts in such a way that small quantities 
of gas are emitting at very high temperatures, thus boost- 
ing the [O in] A4363 line. Shocks are of course one of the 
options (Peimbert et al. 1991, Martin 1996, 1997), as well 
as conductive heating at the interface of an X-ray plasma 
with optically visible gas (Maciejewski et al. 1996). Such 
ideas need to be examined quantitatively, and applied to 
the case of I Zw 18, which we shall attempt in a future 
work. 

What are the consequences of our failure in under- 
standing the energy budget on the abundance determi- 
nations in I Zw 18? It depends on how the electron tem- 
perature is distributed in the ++ zone. As emphasized by 
Peimbert (1967, 1996) over the years (see also Mathis et al. 
1998), the existence of even small zones at very high tem- 
peratures will boost the lines with high excitation thresh- 
old like [O in] A4363, so that the temperature derived from 
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[O in] A4363/5007 will overestimate the average temper- 
ature of the regions emitting [O in] A5007. Consequently, 
the true O /H ratio will be larger than the one derived by 
the standard methods. The C/O ratio, on the other hand, 
will be smaller than derived empirically, because the ultra- 
violet C in] A 1909 line will be extremely enhanced in the 
high temperature regions. Such a possibility was invoked 
by Garnett et al. (1997) to explain the high C/O found in 
I Zw 18 compared to other metal-poor irregular galaxies. 
Presently, however, too little is known both theoretically 
and observationally to estimate quantitatively this effect, 
and it is not excluded that the abundances derived so far 
may be correct within 30%. Obviously, high spatial reso- 
lution mapping of the [O ill] A4363/5007 ratio in I Zw 18 
would be valuable to track the origin of the high [O III] 
A4363 seen in spectra integrated over a surface of about 
10 " 2 . 

Beside demonstrating the existence of a heating prob- 
lem in I Zw 18, our photoionization model analysis led to 
several other results. 

The intensity of the nebular He n A4686 line can be 
reproduced with a detailed photoionization model having 
as an input the stellar radiation field that is consistent the 
observed Wolf-Rayet features in I Zw 18. This confirms the 
results of De Mello et al. (1998) based on simple Case B 
recombination theory. 

By fitting the observed [O m]/[0 II] ratio and the an- 
gular size of the NW component with a model where the 
stellar radiation flux was adjusted to the observed value, 
we were able to show that the H II region is not ionization 
bounded and about half of the ionizing photons are leak- 
ing out of it, sufficient to explain the extended diffuse Ha 
emission observed in I Zw 18. 

While the [O i] emission is not reproduced in simple 
models, it can easily be accounted for by condensations 
or by intervening filaments on the line of sight. There is 
no need to invoke shocks to excite the [O i] line, although 
shocks are probably involved in the creation of the fila- 
ments, as suggested by Dopita (1997) in the context of 
planetary nebulae. 

The intrinsic Ha/H/3 ratio is significantly affected by 
collisional excitation: our photoionization models give a 
value of 3.0, to be compared to the case B recombination 
value of 2.75 used in most observational papers. Conse- 
quently, the reddening is smaller than usually estimated, 
with C(H/3) practically equal to 0. 

Our models can be used to give ionization correction 
factors appropriate for I Zw 18 for more accurate abun- 
dance determinations. However, the largest uncertainty in 
the abundances of C, N, O and Ne ultimately lies in the 
unsolved temperature problem. 

It would be, of course, of great interest to find out 
whether other galaxies share with I Zw 18 this [O III] 
A4363/5007 problem. There are at least two other cases 
which would deserve a more thorough analysis. One is 
the starburst galaxy NGC 7714, for which published pho- 



toionization models (Garcia- Vargas et al. 1997) also give 
[O III] A43G3/5007 smaller than observed. However, it still 
needs to be demonstrated that this problem remains when 
modifying the model assumptions (e.g. the gas density dis- 
tribution, possible heating of the gas by photolelectric ef- 
fect on dust particles etc.). In the case of NGC 2363, 
the photoionization models of Luridiana et al. (1999) that 
were built using the oxygen abundances derived directly 
from the observations yielded a [O in] A4363/5007 ratio 
marginally compatible with the observations. These au- 
thors further argued that, due to the presence of large 
spatial temperature fluctuations, the true gas metallicity 
in this object is higher than derived by empirical methods. 
In such a case, the [O in] A4363/5007 ratio becomes even 
more discrepant. 

It might well be that additional heating sources exist 
in giant H n regions, giving rise to such large temperature 
variations and and enhancing the [O in] A4363 emission. 
As mentioned above, such a scenario needs to be worked 
out quantitatively. 

Further detailed observational and theoretical stud- 
ies of individual objects would be helpful, since we have 
shown that with insufficient observational constraints, 
high [O in] A4363/5007 may actually be produced by 
photoionization models. The effort is worthwhile, since it 
would have implications both on our understanding of the 
energetics of starburst galaxies and on our confidence in 
abundance derivations. 
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Fig. 1. WFPC2 Ha HST image of I Zw 18 aligned on 
NW-SE knots (WF3 image transformed to PC resolution, 
cf. Fig. 1 from De Mello et al. 1998; upper region: NW, 
lower region: SE) The circle (radius 2.5") denotes the area 
used for the measurements of the Ha and F336W flux (cf. 
Sect. 3.1.1, Table 2). The rectangle illustrates the aperture 
used for the spectroscopic observations of Izotov & Thuan 
(1998) centered on the NW knot. The absolute position 
is indicative only. F336W and F555W images indicate a 
small shift (~ 0.5-1 ") between the peak emission and 
the center of the Ha emission (as shown by the circle). A 
shift of 0.5 " was used for the aperture of Izotov & Thuan 
(1998). 
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Fig. 2. Observed [O ill] A4363/5007 ratios in I Zw 18 
and corresponding electron temperature. 1: Lequeux et al. 
(1979); 2: French (1980); 3: Kinman & Davidson (1981); 4: 
Davidson & Kinman (1985); 5: Dufour et al. (1988); 6, 7, 
8: Pagel et al. (1992); 9: Skillman & Kennicutt (1993); 10: 
Stevenson et al. (1993); 11: Martin (1996); 12: Legrand 
et al. (1997). 13: Izotov & Thuan (1998); 14: Vi'lchez & 
Iglesias-Paramo (1998). The error bars are those quoted 
by the observers. 
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Fig. 3. a) g(He+)/Q(H°) for starburst model at Z=0.004 
from De Mello et al. (1998) as a function of age. b) 
Q(He + )/Q(H°) as a function of effective temperature for 
different atmosphere models. SdK97: O dwarf models of 
Schaerer & de Koter (1997) at Z=0.004, SLG92: WR 
atmosphere models of Schmutz et al. (1992). Kurucz 
models have negligible Q(He + ) values, c) Hardness pa- 
rameter Q(He°)/(5(H ) for starburst model as in a. d) 
Q(Hc + )/(5(H°) for same atmosphere models as in b and 
Kurucz (1991) models for [Fc/H]=-1.5 and log<?=5. 
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Fig. 4. Comparison of observed properties of I Zw 18 with 
a series of photoionization models of uniform spheres with 
n= 100 cm~ 3 and a filling factor e=0.01 having, as an 
ionizing source, the starburst models from de Mello et al. 
(1998) at different ages. Left panel: models. Right panel: 
observations, with indication of the tolerance for the mod- 
els. The following quantities are plotted: log F(H/3) + 15 
(open circles), log He II A4686/H/3 + 2 (squares), angular 
radius 9 (crosses), [O in] A4363/5007 x 100 (open trian- 
gles), [S n] A6731/6717 (diamonds), log ([O m]/[0 n]) 
(black triangles), log ([S m]/[S n]) (asterisks) and log 
[O l]/H/3 +3 (plus) as a function of the starburst age. 
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Fig. 5. Same as Fig. 4 for a series of models of uniform 
spheres, with n— 100 cm~ 3 and different filling factors, 
ionized by the radiation field from the model of de Mello 
et al. (1998) at an age of 3.1 Myr. 
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Fig. 6. Same as Fig. 5 for a series of models having a shell 
geometry (6 in =1.5"). 
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Fig. 7. Cumulative surface brightness profiles (integrated 
up to radius r) as a function of r 2 in arcsec. All profiles are 
normalised to r — 2.5". Solid lines: models for spherical 
shells of constant density and various inner radii. Dashed 
line: observed Ha profile from the HST image shown in 
Fig. 1. The observations show a clear deviation from spher- 
ical gas distributions. 
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Table 4. Mean relative ionic abundances and mean ionic temperatures for the best uniform sphere model of Fig. 5 





I 


II 


Ill 


IV 


V 


VI 


VII 


VIII 


IX 


H 


2.585e-03 
16029. 


9.974e-01 
17136. 
















He 


3.963e-03 
16089. 


9.658e-01 
16892. 


3.023e-02 
24967. 














C 


7.775e-05 
15846. 


1.074e-01 
16168. 


7.930e-01 
16947. 


9.525e-02 
19230. 


4.243e-03 
29282. 










N 


2.012e-04 


1.047e-01 


7.823e-01 


1.115e-01 


1.090e-03 


1.863e-04 










15621. 


15968. 


16927. 


19533. 


28682. 


35647. 








O 


4.044e-04 


1.105e-01 


8.750e-01 


1.291e-02 


1.126e-03 


4.627e-05 


6.298e-06 








15594. 


15989. 


17126. 


26194. 


31082. 


36044. 


40715. 






Ne 


4.870e-05 


4.810e-02 


9.367e-01 


1.282e-02 


2.286e-03 


8.458e-05 


5.971e-07 








15763. 


16173. 


17035. 


25411. 


30494. 


35832. 


41693. 






Mg 


6.431e-04 


6.247e-02 


8.998e-01 


2.928e-02 


7.267e-03 


5.240e-04 


1.084e-05 








16023. 


16154. 


16904. 


23329. 


27900. 


32795. 


38679. 






Si 


9.374e-05 


8.091e-02 


7.370e-01 


1.374e-01 


4.409e-02 


5.123e-04 










16170. 


16314. 


16779. 


17814. 


22287. 


30426. 








S 


3.276e-06 


2.150e-02 


5.415e-01 


4.010e-01 


3.550e-02 


4.008e-04 


6.938e-05 








15734. 


15858. 


16314. 


17935. 


21174. 


28651. 


37285. 






CI 


1.374e-05 


4.533e-02 


5.734e-01 


3.636e-01 


1.670e-02 


9.194e-04 


3.075e-05 


9.525c-06 






15835. 


15938. 


16396. 


18119. 


23487. 


29260. 


34968. 


40585. 




Ar 


2.431e-06 


5.553e-03 


6.254e-01 


3.578e-01 


9.725e-03 


1.437e-03 


8.322e-05 


3.004e-06 


7.254e-07 




15517. 


15865. 


16329. 


18257. 


26182. 


29982. 


33949. 


39752. 


43983. 


Fe 


8.416e-06 


1.722e-03 


8.484e-02 


8.665e-01 


3.197e-02 


1.204e-02 


2.750e-03 


2.036e-04 


3.611e-06 




15512. 


15518. 


15860. 


16894. 


22399. 


26289. 


30082. 


34962. 


39424. 



Table 5. Mean relative ionic abundances and mean ionic temperatures for the best spherical shell model of Fig. 6 





I 


II 


Ill 


IV 


V 


VI 


VII 


VIII IX 


H 


2.410e-03 
16315. 


9.976e-01 
16709. 














He 


3.938e-03 
16324. 


9.651e-01 
16658. 


3.092e-02 
18322. 












C 


7.507e-05 


1.109e-01 


8.212e-01 


6.743e-02 


3.274e-04 










16188. 


16392. 


16710. 


17188. 


18455. 








N 


1.532e-04 


1.011e-01 


8.224e-01 


7.626e-02 


1.352e-04 


2.690e-07 








16013. 


16271. 


16714. 


17216. 


18287. 


18819. 






O 


2.960e-04 


1.073e-01 


8.868e-01 


5.670e-03 


2.764e-05 










15992. 


16277. 


16749. 


18397. 


18731. 








Ne 


4.440e-05 


5.038e-02 


9.426e-01 


6.959e-03 


5.394e-05 


0. 


0. 






16110. 


16381. 


16714. 


18208. 


18616. 


0. 


0. 




Mg 


6.280e-04 


6.240e-02 


9.202e-01 


1.652e-02 


2.233e-04 


4.500e-07 


0. 






16285. 


16371. 


16705. 


18116. 


18547. 


18831. 


0. 




Si 


1.011e-04 


7.953e-02 


7.804e-01 


1.244e-01 


1.549e-02 


1.379e-05 








16431. 


16467. 


16676. 


16956. 


17567. 


18010. 






S 


3.123e-06 


2.051e-02 


5.821e-01 


3.819e-01 


1.545e-02 


4.877e-05 


2.998e-08 






16083. 


16189. 


16466. 


17062. 


17737. 


18475. 


19133. 




CI 


1.317e-05 


4.477e-02 


6.190e-01 


3.293e-01 


6.927e-03 


6.193e-05 


8.940e-08 


0. 




16195. 


16247. 


16516. 


17102. 


18047. 


18668. 


19015. 


0. 


Ar 


1.750c-06 


5.340e-03 


6.789e-01 


3.139e-01 


1.878e-03 


2.086e-05 


3.481e-08 






15837. 


16194. 


16499. 


17159. 


18389. 


18745. 


19137. 




Fe 


5.725e-06 


1.128e-03 


7.662e-02 


8.766e-01 


4.300e-02 


2.612e-03 


3.737e-05 


2.083e-08 0. 




15901. 


15917. 


16173. 


16676. 


18213. 


18661. 


18797. 


19168. 0. 



